Edge states of Spin-1/2 Two-Leg Ladder with Four-Spin Ring Exchange 
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A topological insulator and its spin analog as a gapped spin liquid have characteristic low energy 
excitations (edge states) within the gap when the systems have boundaries. This is the bulk-edge 
correspondence, which implies that the edge states themselves characterize the gapped bulk spin 
liquid. Based on the general principle, we analyzed the vector chirality and rung singlet phases of 
the spin- \ ladder with ring exchange by using the edge states and the entanglement entropy. 
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I. INTRODUCTION 



Spontaneous symmetry breaking has been fundamen- 
tal to characterize phases of matter. An example is a 
magnetic phase where a local order parameter is defined 
as a quantum or thermal average of local combination of 
spins. Its ordered phase is not invariant against the sym- 
metry operation which leaves the Hamiltonian invariant. 
In spite of enormous success of this concept, it has been 
also realized that there still exist many important phases 
that cannot be captured by the spontaneous symmetry 
breaking. Strong quantum fluctuation in low dimensions 
prevents such formation of the ordered state and makes it 
possible to realize states without any fundamental sym- 
metry breaking. Such phases form a novel class of matter 
as quantum liquids and spin liquids. Typical examples of 
the quantum liquids are many of integer and fractional 
quantum Hall states. Spin analog of the quantum liq- 
uids is the spin liquid which includes the Haldane integer 
spin chainii 2 ., generic valence bond solid (VBS) states 3 -^, 
and some of the exactly solvable models that also have 
these spin liquid ground states (Kitaev model^, tensor 
category^ and so on). 

As for an excitation of a quantum system, local char- 
acters of the order parameter and the symmetry break- 
ing are fundamental. Gcnerically speaking, one requires 
some mechanisms to have a gapless excitation. Typi- 
cal machineries are the existence of the Fermi surface 
(formed by fermionic quasiparticles) and the Nambu- 
Goldston bosons associated with the broken continuous 
symmetry. An example of the latter is the spin wave 
excitation associated with the Neel order. Gcnerically 
one may construct a gapless excitation by spatial mod- 
ulation of the local order parameter such as the Lieb- 
Shultz-Mattis type spin twists^. In contrast, as for the 
quantum liquids and the spin liquids, it is natural to have 
a finite excitation gap unless one assumes exotic spinon 
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Fermi surfaces. 

To characterize such quantum liquids and the spin liq- 
uids, topological quantities such as the Berry phases and 
the Chern numbers are quite useful^. They do not 
have any symmetry breaking nor local order parame- 
ter. Still, there are many kinds of interesting quantum 
states without characteristic low energy excitation. Such 
non-trivial quantum phases are the topological insulator 
and its spin analog. A generic concept to discuss such 
a phase is the topological order which was first consid- 
ered for the quantum Hall states^. As in the quantum 
Hall state a 11 ' 12 i 13 i 14 ; non-trivial topological insulators do 
have characteristic edge states or impurity states. The 
bulk is gapped and insulating. However the existence of 
boundaries or impurities brings low energy excitations. 
The quantum Hall state of graphene also belongs to this 
topological insulato r 15 i 16 ' 17 . This is the bulk-edge corre- 
spondence where topologically non-trivial bulk guaran- 
tees the existence of localized modes and such low en- 
ergy localized excitations characterize the gapped bulk 
insulator— conversely. 

Not only in the electronic systems, but also, in a quan- 
tum spin system such as the Haldane spin chain, is the 
bulk-edge correspondence realized as the existence of the 
Kennedy triplet for an open chaini^, which was confirmed 
experimentally as welt^. As a novel theoretical tool, the 
entanglement entropy (E.E.) for the gapped topological 
insulators has been quite successful to classify the VBS 
states of the spin chains^^, where this bulk-edge corre- 
spondence plays a fundamental role. 

The spin-i two-leg ladder with four-spin ring exchange 
has rich phase structure due to the frustration 2 ^. For this 
model, the static and dynamical properties have been 
intensively rt u & e &L22&2*&2S2L22&. Although the 
ring exchange model is simple, due to the frustration, 
the ground state is quite involved, namely in the vector- 
chirality (VC) phase. In this paper we consider the VC 
phase and the rung-single (RS) phase of this model to 
identify the ground state properties from the view point 
of the bulk-edge correspondence. 
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FIG. 1: Lattice structure: J is the exchange interaction in the 
leg or rung direction and K is the four-spin cyclic interaction. 
x-(y-) axis is defined as the leg (rung) direction. 
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FIG. 2: Two types of the dimerized plaquette-singlets states. 
Shaded links denote the singlet of the two 5=1/2 spins. 



II. MODEL 



The Hamiltonian is given by 



j 



N/2 N/2 
x — 1 y— 1,2 x — 1 



N/2 

+A"V(P K +P- 1 ). 



x = l 



(1) 



Hereafter, we parametrize the exchange parameter (J, K) 
as (J,K) = (cos 9, sin 6). The four spin cyclic exchange 
consists of the two and four spin exchange interactions 



as' 
P, 



+ 8^,2 • Sx+1,2 + S X) 1 • Sj; + i : 2 + S x ,2 • &x +1,1 

+4(S a;j i • S a; +i i i)(S a ;+i ) i • S-j+i^) 

-4(8,,! • S x+li2 ){S Xt2 • Sa+i.1) + 1/4. (2) 



Here, S x , v denotes a spin-i operator at site (x,y) (see 
Fig. [I}. We assume the system-size N to be even. 



III. EDGE STATES 

The Berry phases as topological order parameters clas- 
sify the rung-singlet and the VC and RS phases of the 
model 2 ^. 

In the VC phase, Hamiltonian |T]) is adiabatically 
connected to a decoupled vector-chiral model Hdvc = 
IL=ovon( S z,i x S x,2) ■ ( s z+i,i x Sx+1,2), whose ground 
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FIG. 3: Two types of open boundary conditions (OBCs). 

state is a direct product of plaquette-singlets \ip a ) [see 
Fig. [Ha)] , which is defined as 

|V0 = [(1,1), (2, 2)]® [(1,2), (2,1)] 

®[(3,1),(4,2)]®[(3,2),(4,1)]®-" (3) 

Here we have introduced a singlet at two sites a and /?, 
[a,P] = [(a x ,a y ),{f3 x ,Py)] as 



[S,0\ 



1 
71 



(|T>s|i>£-|T>0| i)a) 



(4) 



By taking into account of the translational invariance, 
we consider the linear combination of \ip a ) and as 
shown in Fig. HJb). State |-06) is defined as 



\A) = [(2,1), (3, 2)]® [(2, 2), (3,1)] 

®[(4,1),(5,2)]® [(4,2),(5,1)]® 



(5) 



Notice that the two dimerized plaquette singlet states 
\t/) a ) and are not orthogonal to each other, but their 
overlap is exponentially small for large N. For N/2 being 
even, the overlap is obtained as (V'alV'fc) = 2 2 ~ w / 2 . As 
described later, the state \ip s ) oc \ip a ) + \ipi) can be a 
good trial state to understand the edge states and the 
entanglement entropy for the vector chirality stated. 

In the VC phase, the bulk itself has a finite gap in 
the thermodynamic limi t 21 i 24 . We introduce two types 
of boundaries — (a) diagonal-edge and (b) vertical-edge 
(see Fig. [Hp 2 -. We diagonalizcd the Hamiltonian numer- 
ical by the Lanczos method. In Fig. 2] we show the size 
dependence on the energy gap of the total S z = sec- 
tor for 6 = 47r/5, where the bulk spin gap is relatively 
large. The system with the vertical-edges has almost the 
same excitation energy as that of the periodic one. In 
contrast, there exist additional low energy excited states 
with an exponentially small energy gap when the sys- 
tem has the diagonal-edges. The appearance of such a 
localized mode indicates the feature of the bulk-edge cor- 
respondence. In fact, this mode is a triplet excitation and 
can be interpreted as the Kennedy triplet (see Fig. [5]) - 
triplet excitation between the effective boundary spins 
at both sidesi^. We confirmed that this mode is really a 
triplet excitation from the diagonalization on the differ- 
ent S z sector. This mode is observed also in the different 
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FIG. 4: System-size dependence of low excitation energy for 
9 = 4-7r/5 with the two OBCs compared with the energy gap 
under the periodic boundary condition (PBC). 
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FIG. 5: Schematic picture of the low energy excitations. 



(J,K)'s in VC phased. Note that the Kennedy triplet 
in this model was also discussed in a different context 
before^. 

When we introduce the boundary for the trial state 
\ip s ), i- e - ; the linear combination of the two dimerized 
plaquette singlets, the isolated spins appear near the 
boundaries (see Fig. [5]). In the system with diagonal- 
edges they appear at each boundary while in the system 
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FIG. 6: Linear combination of the dimerized plaquette- 
singlets \ip a ) with OBCs. 



with vertical-edges they appear as a pair at one side of the 
edge. Although a pair of localized spins behaves freely in 
the decoupled vector chiral model 7idcv, the pair of the 
spins couples each other in the original model 7i C yc with 
vertical-edges. Therefore, the Kennedy triplet excitation 
does not appear in the system with the vertical-edges. 
Thus the trial state \ip a ) gives consistent understanding 
for the low energy spectra of the VC phase. 

Next we consider the RS phase. In Fig. [7j we show the 
size dependence on the energy gap of the total S z = 
sector for 6 = —tt/5. Combining with the calculation on 
different S z sector we obtain the Kennedy triplet mode 
as in the VC phase in the case with diagonal-edge. This 
is consistent with the naive picture of the rung singlet, 
which is obtained by the Berry phase (see Fig. [SJ . This 
mode is confirmed also in the different (J, K)'s in the RS 
phase. 



IV. ENTANGLEMENT ENTROPY 

In this section, let us consider entanglement entropy for 
the VC state. Define the density matrix of the state 
as p = \ip){ip\ and divide the system into two subsystems 
A and B [see Figs.[9ta) and[9^b)]. Then the entanglement 
entropy is defined as^ 



E.E. 



-(log /3a) 



-Tr [p A log p A ] 



Here the reduced density matrix pa is given as pa = 
TtbP- It represents how much state \tp) is entangled be- 
tween subsystems A and B. 

Let us assume state as a unique ground state un- 
der the PBC in the VC state. We take subsystem A as 
the subsystem with the diagonal (vertical) edges with N' 
spins by the reduction from N — 20 spin system with 
the PBC. Similarly as for the energy gap we consider 
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FIG. 7: System-size dependence of low excitation energy for 
9 — — tt/5 with the two OBCs compared with the energy gap 
under the periodic boundary condition (PBC). The inset is 
the extended figure for the first excitation state energy. 
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FIG. 8: Rung-singlets with OBCs. 



two types of open boundaries, we calculate the entangle- 
ment entropy numerically for the cases with vertical-edge 
and diagonal-edge as shown in Figs. [S^a) and[9jb). Fig- 
ure Mc) shows the N' dependence of the entanglement 
entropy. In both cases, the obtained E.E. contains a con- 
tribution around 3 log 2. The contribution 3 log 2 can be 
understood by the trial state \tf> s ) (see the Appendix). 

In the naive picture, in the RS phase the isolated single 
on the rung is a good trial state. Therefore we expect 
that E.E.= 2 log 2 for the diagonal edge and that E.E. 
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FIG. 9: Two types of the reduction to the subsystem A with 
(a) diagonal-edge and (b) vertical-edge, (c) Entanglement 
entropy of subsystem for 9 = 4ir/5 with N' spins by reduction 
from N = 20 spin system. 



= for the vertical edge. We have qualitative agreement 
that the E.E. for the diagonal edge is larger than that for 
the vertical edge, although there is strong ^-dependence 
and large deviation from the 2 log 2. This is due to the 
relatively large overlap on the leg-direction. 



V. SUMMARY 

In summary, it has been shown that in the vector chi- 
rality state the dimerized plaquette singlet state \ip a ) can 
be a good trial state to understand the numerical results 
for the topological properties — the edge states and the 
entanglement entropy. The entanglement entropy of \ip s ) 
has been obtained as 3 log 2 for both types of the reduced 
systems while the appearance of the edge states depends 
on the type of boundaries. These boundary- dependent 
low energy excitations as the generic edge states charac- 
terize the vector chirality phase. This boundary depen- 
dent edge state is also observed in the rung-singlet phase. 
Such localized modes in the boundary arc expected to be 
observed experimentally through the impurity or surface 
effect. 
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FIG. 10: Two types of the tensor products of the two dimer- 
ized states in Eq. (|A6|I with diagonal-edges for even N' /2 case 
(N' — 12), which corresponds to Fig[6](a). 



We introduce the dimerized states \D) and \D aa r) on the 
chain with length N' /2, 



\D) = [l,2](g)[3,4]®-..(g)[JV72-l,JV72]^A4) 
\D) aa . = [2,3]®--- 

®[N'/2- 2,N'/2- 1] ® \a') N , /2 . (A5) 
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APPENDIX A: ENTANGLEMENT ENTROPY 
FOR THE TRIAL STATE 



Changing site indexes to decouple the ladder into two 
chains, we have equivalences l-A^o-/) = \D) ® \D) aa i and 
\B aa >) = \D) aa : ® \D) (see Fig. ITU]). Here we have intro- 
duced the normalized state \D) oc \D) — (D^i\D)\D-\{) — 
(D^\D)\D^) by the Gram-Schmidt orthogonalization 
method. Then we have the following relation: 



For example, we calculate the entanglement entropy 
of \ip s ) in the diagonal-edge case for even N/2 and N'/2 
cases. The reduced density matrix pa can be obtained 
as (see Fig. ITU|) . 



PA 



\B„ 



[\A n ){B n \-\A n ){B^\ 



( A ~ B)1A1) 



where 
of the 



\B„ 



m = 2 + 2^ N ' 2 and m! = 2 1 -( JV " W ')/ 2 . States 
and \B aa i) are represented as the tensor product 
dimerized states in the zigzag chains (see Fig. [10]) , 

) = [(l,l),(2,2)]®|a) (2il) 

®[(3,1),(4,2)]®[(3,2),(4,1)]®--. 
®[(Y72-1,1),(Y72,2)] 
®[(iV72 - 1, 2), (N'/2, 1)] ® \a') {N , /2+1>2) , 

(A2) 

) = |a) (2jl) ®[(2,l),(3,2)]®[(2 ! 2),(3,l)]®--- 
®[(Y72-2,1),(Y72-1,2)] 
®[(N'/2 - 2, 2), (N'/2 - 1, 1)] ® \a') {N , /2 , 2) 
®[(N'/2,l),(N'/2+ 1,2)]. (A3) 



p A ~ - ^ III- )® \Daa') )[(D\® (D aa , 
<r,<r'=U 



(7,0 

+ {\D 



(A6) 



The approximation holds up to the order 0(2 N I 2 ) + 
0{2~( N ~ N )/ 2 ) and the eight summands correspond to 
the states shown in Fig. ]6](a). Thus, the entanglement 
entropy is obtained as E.E. ~ —8 x glogg = 3 log 2. 
In a similar manner we can calculate the entanglement 
entropy in the case of the other geometry. The entangle- 
ment entropy counts the degrees of freedom of the spins 
around the edges, although the Kennedy triplet does not 
appear in the system with vertical-edge due to the short 
range residual interaction between the localized effective 
spins. 
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